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EXECUTIVE  SUMMARY 


OBJECTIVE 

One  of  the  High  Temperature  Superconductor  (HTS)  bi- Superconducting  Quantum  Interference  De¬ 
vice  (SQUID)  Array  Fractional  Wavelength  Baseline  Determination  and  Design  Study  team’s  goal  was 
to  determine  whether  the  raw  simulated  output  from  an  array  of  bi-superconducting  quantum  interference 
devices  (bi-SQUIDs)  could  perform  direction-finding  (DF)  tasks  such  as  the  angle  of  arrival  (AoA),  eleva¬ 
tion,  and  polarization.  First,  we  applied  the  multiple  signal  classification  (MUSIC)  algorithm  on  linearly 
polarized  signals.  We  included  multiple  signals  in  the  output  both  of  the  same  frequency  and  different  fre¬ 
quencies.  Next,  we  explored  a  modified  MUSIC  algorithm  called  dimensionality  reduction  MUSIC  (DR- 
MUSIC),  which  can  handle  elliptical  or  circular  polarizations.  Finally,  we  explored  the  six-element  Poynt- 
ing  vector  method. 

RESULTS 

We  determined  that  the  MUSIC  algorithm  is  able  to  determine  the  AoA  from  the  simulated  SQUID 
data  for  linearly  polarized  signals.  The  MUSIC  algorithm  could  accurately  find  multiple  signals  of  differ¬ 
ent  frequency,  and  if  there  were  k  +  1  sensors,  it  could  find  k  signals  of  the  same  frequency.  Unfortunately, 
if  signal  polarizations  arc  circular  or  elliptic,  desired  phase  differences  (due  to  the  antenna  separations)  arc 
changed  by  the  signal  polarizations  and  the  MUSIC  algorithm  cannot  find  the  direction-of-arrival  (DoA) 
angles.  We  explored  some  modified  algorithms,  including  DR-MUSIC,  which  performed  well  with  both 
smaller  baselines  and  signal  polarizations. 

The  Poynting  vector  algorithm  performed  well  for  any  polarization.  The  disadvantage  of  this  method 
is  that  it  requires  electric  field  antenna  measurements  to  capture  all  six  components.  As  a  result,  the  SQUID- 
based  sensor  size  benefit  is  reduced  and  applications  are  restricted.  Outside  frequencies  of  interest,  the 
Poynting  vector  error  increased  when  the  separation  between  arrays  was  the  same  magnitude  as  the  signal 
wavelength. 

RECOMMENDATIONS 

We  determined  it  is  possible  to  apply  DF  when  using  bi-SQUID  array  sensors.  Future  tests  should 
scale  up  the  array  used  and  determine  the  distance  limits  between  arrays.  We  are  also  exploring  other  con¬ 
figurations  that  use  both  phase  difference  and  Poynting  vector  algorithms  for  all  DF  tasks  with  any  polar¬ 
ization  without  the  need  to  use  electric  field  sensors.  Finally,  experiments  in  this  report  were  conducted 
using  computer  simulation  of  a  bi-SQUID  array  sensor;  we  will  need  to  verify  all  algorithms  with  experi¬ 
mental  data. 
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1.  INTRODUCTION 


Superconducting  quantum  interference  devices  (SQUIDs)  arc  magnetometers  that  use  macroscopic 
quantum  effects  to  achieve  unprecedented  sensitivity  to  surrounding  magnetic  fields.  A  direct  current  (DC) 
SQUID  has  two  Josephson  junctions  in  a  ring  of  superconducting  material.  Adding  a  third  junction  bi¬ 
secting  the  superconducting  loop  increases  the  lineality  of  the  voltage  response  output  [1],  This  device  is 
called  a  DC  bi-SQUID. 

Through  recent  developments  [2,  3],  the  idea  of  turning  a  SQUID  (or  bi-SQUID)  or  an  array  of  SQUIDs 
(or  bi-S  QUIDs)  into  a  working  antenna  is  becoming  a  reality.  One  step  to  accomplishing  this  goal  is  to 
adopt  or  develop  a  method  that  would  allow  a  SQUID  (or  bi-SQUID)  array  sensor  to  perform  direction 
finding  (DF)  (i.e.,  locate  the  source  of  a  signal).  The  simplest  way  to  achieve  this  performance  goal  is  to 
spatially  separate  multiple  SQUID  (or  bi-SQUID)  arrays  to  exploit  natural  phase  differences  induced  by 
separations.  Other  methods  involve  distinctly  polarized  arrays  that  sense  the  Cartesian  components  of  the 
electromagnetic  waves. 

In  this  report,  we  apply  apply  multiple  signal  classification  (MUSIC)  [4],  dimensionality  reduction 
MUSIC  (DR-MUSIC)  [5],  and  Poynting  vector  [6]  algorithms  to  a  small  simulated  array  of  bi-SQUID  ar¬ 
ray  sensors  and  record  their  performance  against  a  variety  of  different  arrival  angles  and  polarizations.  To 
maximize  the  utility  that  a  bi-SQUID  array  may  offer,  the  project  team  focused  on  small  and  linear  arrays 
consisting  of  two  to  three  antennas  so  that  they  arc  installed  on  different  platforms.  We  are  exploring  other 
potential  methods  [7-11]. 
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2.  BACKGROUND:  ARRAY  SIMULATION 


To  simulate  data  that  is  fed  into  the  DF  algorithms,  we  use  the  schematic  of  the  series  coupled  array  of 
DC  bi-SQUIDs  in  Figure  1  to  derive  a  system  of  equations  that  model  the  Josephson  junction  phase  differ¬ 
ences.  In  the  circuit,  (4,  4,fc,  %2 ,k,  h,k,  H,k,  4,fc)  represent  the  normalized  currents,  (pi,k,  P2,k,  P3 ,&))  are 
the  phases  across  the  Josephson  junctions,  ( L\a L]b.k,  L‘>a.k,  L-2b,k ,  ®3a,/b  are  the  parameters  re¬ 
lated  to  the  inductance  values,  and  k  =  1, N,  and  xe  are  the  points  in  the  array  where  the  contributions 
from  external  field  and  input  signals  arc  included. 


To  derive  equations,  we  find  current  equations  for  the  schematic  diagram  using  Kirchhoff’s  current 
law  and  use  the  resistively  shunted  junction  (RSJ)  model  of  the  over-damped  Josephson  junction  with  a 
lumped  circuit  representation.  Both  equations  for  bi-SQUID  arrays  and  their  derivations  are  found  in  [12, 
13].  Some  fabrication  considerations  arc  addressed  in  [14-16].  The  system  of  equations  is 


,  T  ,  ,dp l,fc  T  dtp2,k  T  d^3:k 

\-^l,k  +  1-J2a,k)  -LJ2b,k  1  -^1  ,k  1 

dr  dr  dr 

—  b,k^B  H-  ^e,k  ^P2,k  ^1, k  (-^1,/c  “1“  -^2a,/c)^cl,fc  sin  ^  +  L\,kijc?>,k  sin  ^Ps,k 

T  •  •  \  , 

“1“  IJ2b,k?'c2,k  9^2, /c  /  J  "Tg  9^2,2  2i7TXeCli) 

i^k  di  ai 

T  dip  t;fc  .  ^dp2,k  T  d(f3tk 

L2d— - (Li  k  +  ^2b,k)—i - -Fi, A: — - - 

ar  dr  dr 

—  Lla,k^B  T  Pe,k  ^P2,k  ^Pl,k  T  ( Ll,k  T  dJ2b,k)ic2,k  sill  <P2,k  4~  L\ ,k^c3,k  sin  P3,k 

T  •  •  \ ^  /  C\  \ 

IJ2a,k'l/cl,k  Sill  +  y  ^  (ty?l,2  ^2 ,i  27 TXeCli) 

i^k  di  ai 

T  difi’k  dp2,k  T  dif3yk 

X'2a,k  j  lJ2b,k  r  -^3 ,k 

dr  dr  dr 

=  P2 ,k  -  <Pl,k  —  <F3 ,k  —  L2a, kid, k  sin  Pltk  +  L  2b, kic2,k  SHI  C p2,k  ~  ^3,fc*c3,fc  sin  ip3tk , 


(1) 


where  ipi}k,  P2,k,  and  993  k  are  phases  across  the  Josephson  junction  in  each  bi-SQUID,  k  =  1, ...,  N, 

N  is  the  number  of  bi-SQUIDs  in  the  array,  M  is  the  coupling  parameter,  and  L\k  =  Tia.fc  +  Lib,k, 

Lz,k  =  d,3a,k  +  T 3 Dots  denote  the  time  differentiation  with  normalized  time  r  =  uct,  where  t  is 
time,  ujc  =  2eI()Rr\i / h,  and  the  normalizing  current  of  the  Josephson  junctions  is  Jo-  The  parameter  II ^ 
in  uic  is  the  normal  state  resistance  of  the  Josephson  junctions,  e  is  the  charge  of  an  electron,  and  h  is  the 
reduced  Planck  constant.  The  parameter  4  is  the  normalized  bias  current  and  tpe  k  =  2irxeak,  where  ak  is 
the  normalized  bi-SQUID  area.  We  use  the  approximate  assumption  that  ak  =  L  \„j-  +  L  \  b , /,■ + L 2 a , / • + L 2 h . /,- ■ 

To  model  the  array,  the  differential  equations  in  Equation  (1)  arc  integrated  in  MATLAB®  and  the 
array  voltage  over  time  V (t)  =  YliLi  ( js  plotted.  See  Figure  2  for  an  example  of  V (t). 
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Raw  Output 


Moving  Average  Output 


Figure  2.  A  bi-SQUID  array  V(t)  with  multiple  signals  of  different  strengths  and  polarizations.  The 
top  plot  is  raw  data  and  the  bottom  is  the  data  after  passing  through  a  filter  that  removes  the 
terahertz  bi-SQUID  response. 


Multiple  sensors  are  necessary  for  the  feasibility  of  any  DoA  estimation.  The  bi-SQUID  simulation 
program  is  only  equipped  to  emulate  a  single  bi-SQUID  array  sensor.  To  compensate,  we  assume  that  we 
have  two  or  more  identical  bi-SQUID  array  sensors  that  are  located  at  distances  we  select.  Inputting  a  sig¬ 
nal  into  the  bi-SQUID  program  results  in  an  estimate  of  the  voltage  response  that  a  true  bi-SQUID  array 
would  create  in  the  presence  of  the  signal.  To  replicate  the  output  of  additional  sensors  that  were  physi¬ 
cally  displaced,  we  input  the  same  signal  but  doped  with  a  predetermined  phase  difference  that  is  specific 
to  the  AoA  we  elected  to  simulate. 

For  K  signals  (si(t), . . . ,  s a'(Q)  impinging  on  a  uniform  linear  array,  the  signal  received  by  the  mth 
antenna  is  given  as 


cfcSfc(f)exp 


/  j2(m  —  1)7tA 

V  ^ 


COS 


(2) 


where  A  is  the  length  of  separation  between  antennas  (baseline),  Xk  is  the  wavelength  of  the  A:lh  signal, 
and  0  <  (j)k  <  2-7T  is  the  arrival  azimuth  of  the  klh  signal.  The  parameter  c/,:  represents  a  complex  scalar 
that  is  derived  from  the  relationship  between  signal  parameters  and  antenna  polarizations.  When  a  bi- 
SQUID  array  simulator  does  not  accept  a  complex- valued  input,  Equation  (2)  is  replaced  with 


Zk&k 


+  dk  + 


j2(m  —  1)7tA 
^ k 


(3) 


where  Zk  =  |c|  and  dk  =  arg(c/,,).  We  can  replace  Equation  (2)  with  Equation  (3)  because  the  exponential 
and  complex  terms  in  Equation  (2)  represent  phase  differences  that  we  can  easily  transfer  into  the  real¬ 
valued  signal  itself.  If  our  signal  takes  the  form  of  Equation  (3),  a  Hilbert  transform  is  applied  to  the  re¬ 
sulting  output  of  the  bi-SQUID  array  simulation.  The  Hilbert  transform  allows  us  to  complete  the  analytic 
signal  necessary  for  use  in  DF  algorithms. 
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3.  APPLYING  THE  MULTIPLE  SIGNAL  CLASSIFICATION  ALGORITHM 


The  first  algorithm  we  explore  is  the  MUSIC  algorithm  [4]  with  a  uniform  linear  array  of  bi-SQUID 
array  sensors. 

3.1  SUMMARY  OF  MUSIC 

Integer  K  number  of  polarized  narrowband  electromagnetic  waves  traveling  through  the  far  field  through 
a  homogeneous  isotropic  medium  impinge  on  an  array  consisting  of  integer  L  number  of  uniformly  polar¬ 
ized  antennas  uniformly  spaced  at  no  more  than  a  half-wavelength  apart  from  each  other.  Data  collected  by 
L  antennas  are  represented  by 


m 


I< 

^afcs  k(t)  +n  (t) 
k= i 


As(f)  +  n(t), 


(4) 


where 


a  Ofc) 


1 

gJ’Mfc 

^3  l)Mfc 


is  the  steering  vector  for  the  kth  signal,  n(t)  is  the  zero-mean  additive  complex  Gaussian  noise  and  ///,.  = 
2ir  A  cos((/>fc) / \k.  For  the  matrix  notation,  A  =  [ai , . . . ,  a/c] .  The  data  matrix  X  is  an  L  x  N  matrix  where 
N  is  the  number  of  data  samples  recorded  in  a  given  interval  of  time. 


The  covariance  matrix  of  the  data  X(f)  is  given  by 


It,  =  E[XXh]  =  ARsAh  +  cr2I, 


where  (-)H  denotes  complex  conjugate  transpose,  a2  is  the  white  noise  power,  and  Rs  =  E[sH(t)s(f)]  is  the 
source  covariance  matrix.  We  can  easily  estimate  R;c  with 

R  =  4XX» 


MUSIC  begins  by  executing  an  eigendecomposition  on  the  matrix  R.  The  resulting  eigenvectors  de¬ 
compose  the  data  correlation  matrix  into  a  if -dimensional  signal  subspace  and  an  (L  —  K) -dimensional 
noise  subspace.  Let  En  represent  the  L  x  (L  —  K)  matrix  of  noise  eigenvectors  associated  with  the  ( L  —  K ) 
smallest  eigenvalues  of  R.  The  MUSIC  spectrum  is  then  given  as 


Pmu(^,  4>k) 


aH(fo)EwEffa(fo) 

(4>k)aH  (4>k) 


The  peaks  of  this  spectrum  will  correspond  to  the  desired  azimuth  angles  to  be  found  (Figure  3).  The  left 
plot  is  the  raw  V (t)  data  with  two  signals  of  300  MHz  and  3  GHz.  The  right  plot  is  the  AoA  for  the  300- 
MHz  signal  as  determined  by  the  MUSIC  algorithm. 

For  a  distinctly  polarized  or  uniform  rectangular  array  of  antennas,  the  corresponding  steering  vec¬ 
tor  a  (Ok,  (pk)  allows  MUSIC  to  perform  a  two-dimensional  (2-D)  search  so  that  elevation  can  be  found  in 
addition  to  the  azimuth.  In  the  case  where  only  azimuth  is  considered,  elevation  is  held  constant  at  7t/2 
radians. 
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Voltage  Response  of  SQUID  simulation  with  300  MHz  and  3  GHz  signals 

i  ?  1  ,  '  I  ’ 


5  ‘  - Original  Data 

Lowpass  Filtered  Data 


2  4  6  8  10  12 

Sample  x104 


Figure  3.  Left:  bi-SQUID  array  V(t)  with  two  signals  of  300  MHz  and  3  GHz.  Right:  the  AoA  for  the 
300-MHz  signal  as  determined  by  the  MUSIC  algorithm. 


3.2  MULTIPLE  SIGNALS,  DIFFERENT  FREQUENCIES 

For  conventional  antenna  arrays  with  a  limited  bandwidth,  there  are  a  limited  number  of  frequencies 
detected.  In  contrast  to  the  limited  bandwidth  of  a  conventional  antenna,  the  bandwidth  of  a  bi-SQUID 
array  sensors  is  extended,  which  leads  to  an  increase  in  active  signals.  To  compensate  for  these  increases, 
two  to  three  bi-SQUID  array  sensors  coupled  with  digital  filters  are  used  to  maintain  their  extremely  small 
profile  while  preserving  acceptable  performance  in  the  presence  of  many  impinging  signals  of  differing 
frequencies.  Digital  filters  offer  the  ability  to  block  unwanted  frequency  ranges  so  that  MUSIC  can  find  the 
DoA  angles  selectively. 

For  the  first  experiment  (Table  1),  we  test  the  MUSIC  algorithm  against  a  bi-SQUID  array  simulation 
of  two  uniformly  polarized  sensors  that  have  received  one  signal.  We  test  limits  by  steadily  decreasing  the 
baseline  until  failure.  The  baselines  listed  are  in  units  of  wavelengths.  Each  test  used  21'  data  samples. 

Table  1.  Azimuth  estimate  using  MUSIC. 


Baseline 

</> 

^est 

Error 

1/2 

145 

145 

0 

1/2 

10 

9 

1 

1/10 

45 

46 

1 

1/10 

120 

119 

1 

1/25 

10 

20 

10 

1/25 

80 

80 

0 

1/25 

165 

156 

9 

1/50 

30 

46 

16 

1/50 

30 

122 

92 

For  the  second  experiment,  we  introduce  one  3-GHz  signal  and  one  300-MHz  signal  into  the  sim¬ 
ulation  and  try  to  find  the  angles  with  MUSIC.  This  experiment  is  for  three  uniformly  spaced  identical 
bi-SQUID  array  sensors.  While  searching  for  both  signals,  it  is  possible  to  tune  the  steering  vector  (i.e., 
change  the  phase  difference  to  match  a  particular  wavelength)  in  our  search  to  focus  on  one  signal  at  a 
time.  In  place  of  this  approach,  we  use  filters  to  block  the  unwanted  signal  and  tell  MUSIC  that  only  one 
signal  is  present.  This  approach  yields  better  results  (Table  2).  Two  baselines  listed  in  each  row  correspond 
to  the  3-GHz  and  300-MHz  signals,  respectively. 
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Table  2.  Two  signals  found  separately  using  filters. 


Baseline 

0  i 

02 

01  est 

02  est 

Error 

1/2,  1/20 

15 

135 

15 

135 

0,  0 

20 

60 

20 

60 

0,  0 

1/5,  1/50 

20 

150 

19 

161 

1,  11 

45 

90 

46 

91 

1,  1 

1/10,  1/100 

110 

160 

110 

180 

0,  20 

20 

55 

23 

87 

3,  64 

3.3  MULTIPLE  SIGNALS,  SINGLE  FREQUENCY 

The  MUSIC  algorithm  has  the  option  of  finding  multiple  signals  of  the  same  frequency  simultane¬ 
ously,  in  addition  to  its  accuracy  and  easy  implementation.  However,  in  the  presence  of  k  signals,  MUSIC 
requires  k  +  I  antennas  for  reliable  DF.  Hence,  for  an  array  with  three  elements,  a  maximum  of  only  two 
signals  with  the  same  carrier  frequency  can  be  located.  For  situations  where  this  constraint  is  too  restric¬ 
tive,  an  array  with  more  sensors  is  preferred.  Table  3  show  results  from  simulations  with  two  incoherent 
signals  impinging  on  a  three-element  uniformly  polarized  array  with  identical  carrier  frequencies.  The 
MUSIC  algorithms’  ability  to  find  signals  simultaneously  is  used  to  find  the  angles  since  is  not  possible  to 
focus  on  one  signal  at  a  time. 


Table  3.  Two  signals  found  simultaneously. 


Baseline 

01 

02 

01  est 

02  est 

Error 

1/2 

110 

160 

109 

158 

1,  2 

40 

140 

35 

134 

5,  6 

15 

115 

8 

114 

7,  1 

1/10 

80 

110 

81 

113 

1,  3 

90 

100 

90 

99 

0,  1 

1/25 

10 

170 

35 

35 

25,  135 

45 

135 

59 

59 

14,  76 

3.4  POLARIZED  SIGNALS 

Unfortunately,  if  signal  polarizations  arc  circular  or  elliptic,  the  desired  phase  differences  invoked  by 
the  antenna  separations  arc  polluted  by  the  independent  phases  that  result  from  the  interactions  of  the  po¬ 
larization  and  the  signal  phase  differences.  The  MUSIC  algorithm  inconsistently  finds  DoA  angles  in  this 
situation  unless  a  modified  technique  is  used  [5,  9,  11].  While  using  an  ordinary  form  of  the  MUSIC  algo¬ 
rithm,  we  restrict  our  signals  exclusively  to  linear  polarization. 
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4.  DIMENSIONALITY  REDUCTION  MULTIPLE 
SIGNAL  CLASSIFICATION 


While  the  previous  section  was  useful  for  signals  with  linear  polarizations,  they  do  not  perform  well 
with  elliptical  or  circular  polarizations.  To  achieve  accuracy  when  working  with  signals  possessing  varied 
polarizational  phase  differences,  we  can  use  a  technique  called  DR-MUSIC  [5].  The  DR-MUSIC  algo¬ 
rithm  is  designed  to  work  with  a  six-element  colocated  vector  sensor.  The  components  of  the  electromag¬ 
netic  wave  of  an  incoming  signal  is  picked  up  by  the  vector  sensor  and  is  expressed  as 


&k(fik i  Tfc)  Vk) 


^ kx 

cos  9k  cos  4>k 

-  sin  4>k 

€-ky 

cos  6k  sin  4>k 

COS0k 

& kz 

—  sin  9  k 

0 

hkx 

-  sin  (j)k 

—  cos  9k  cos  4>k 

hky 

COS  0k 

—  cos  9k  sin  0k 

hkz 

0 

sin  9  k 

(5) 


where  —  it  <  ijk  <  tt  is  the  polarization  phase  difference  of  the  A;lh  signal.  Denote  the  larger  matrix  of 
Equation  (5)  as  f l(9k,  0k)-  The  matrix  of  data  X  is  found  in  the  same  way  as  in  Equation  (4)  but  with  new 
steering  vectors  shown  in  Equation  (5).  Eigendecomposition  is  performed  on  the  data  correlation  matrix 
and  we  now  let  Un  represent  the  noise  eigenvectors  of  the  data.  We  define  a  new  matrix: 


The  DR-MUSIC  algorithm  states  that  the  minimum  eigenvalue  of  B  is  smallest  precisely  when  B  is  evalu¬ 
ated  at  the  correct  DoAl  angles.  Thus,  the  DR-MUSIC  spatial  spectrum  is 


Pdr($Ao  0k) 


-1 


min  Amin(B(0fc,</>fc)) 

8k  i&k 


Table  4  uses  the  modified  version  of  DR-MUSIC  to  find  the  arrival  angles  of  one  signal  impinging  on 
three  orthogonally  oriented  non-colocated  bi-SQUID  array  sensors.  The  signal  is  allowed  to  have  circular 
or  elliptical  polarizations;  therefore,  parameters  are  varied  with  each  trial. 

Table  4.  Non-colocated  DR-MUSIC. 


Baseline 

7 

V 

0 

8 

^est 

^est 

Error 

1/2 

7t/3 

0 

145 

45 

144 

45 

1,  o 

7t/6 

-7t/4 

170 

25 

170 

25 

0,  0 

1/10 

7t/6 

-7t/4 

170 

2 

170 

4 

0,2 

1 

-1 

90 

5 

93 

3.5 

3,  1.5 

1/25 

2tt/5 

7t/6 

10 

90 

9 

90 

1,  o 

2tt/5 

7t/6 

100 

5 

108 

3.5 

8,  1.5 

1/50 

7t/6 

-7t/4 

170 

25 

169 

24.5 

1,  0.5 

2tt/5 

7t/6 

100 

5 

106 

3 

6,  2 

1/100 

tt/2.5 

-37t/4 

10 

20 

8 

18 

2,  2 

7t/5 

7t/4 

12 

87 

12 

86.5 

0,  0.5 

1/500 

7t/5 

7t/3 

80 

80 

80 

79 

0,  1 

7t/3 

7t/2 

145 

45 

146 

45 

1,  o 
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5.  POYNTING  VECTOR 


We  explored  alternate  ways  to  achieve  the  AoA,  elevation,  and  polarization  information  [6-8,  10].  In 
this  section,  we  explore  the  Poynting  vector  method. 

5.1  SUMMARY  OF  THE  POYNTING  VECTOR 

The  Poynting  vector  method  is  the  most  robust  way  to  match  a  small  baseline  bi-SQUID  array  ar¬ 
rangement  with  relatively  long  wavelengths.  The  computations  and  processing  are  simple  and  the  error 
lowers  at  longer  wavelengths.  One  disadvantage  is  that  a  tripole  antenna  is  needed  to  measure  electric 
fields;  therefore,  you  lose  the  advantage  of  using  bi-SQUID  arrays.  This  method  would  also  need  to  dis¬ 
tinguish  multiple  signals  and  handle  multipath  signals.  The  Poynting  vector  method  takes  the  voltage  out¬ 
put  from  each  sensor  to  get  the  field  component  in  that  direction.  Then,  a  cross  product  is  taken  to  find  the 
Poynting  vector  k  (Figure  4). 


Figure  4.  Sample  calculation  of  the  Poynting  vector  from  bi-SQUID  array  data. 

The  direction  of  the  source  is  simply  the  negative  of  the  cross  product: 

k  =  E  x  B. 

One  can  also  do  2-D  DF  with  only  three  sensors:  (Bx,  By,  Ez).  However,  this  only  works  with  vertically 
polarized  sources,  that  is,  when  Bz  =  0.  The  expressions  kx  =  —ByEz  and  ky  =  EZBX  are  used  to  find 
the  x  and  y  components  of  the  Poynting  vector. 

Three  separate  bi-SQUID  array  simulations  are  used  to  generate  magnetic  field  vector  sensor  data, 
each  with  a  specific  input  flux  and  phase  difference  due  to  separation.  A  simple  model  of  electric  fields 
with  some  noise  is  added  to  compute  the  Poynting  vector  and  to  give  the  angle  error  and  vector  quiver 
plots.  The  elevation  and  azimuth  of  the  Poynting  vector  is  the  direction  the  signal  is  traveling,  which  is 
the  negative  of  the  vector  pointing  toward  the  source. 

Equations  for  both  the  x  and  y  components  of  the  magnetic  field  inside  of  a  polarization  ellipse  per¬ 
pendicular  to  the  Poynting  vector  when  pointing  in  the  z  direction  is 

cos  Ob  cos(tuf) 
sin  Ob  cos  (cut  +  ay)  . 

0 
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B  =  B 


Here,  uj  is  the  angular  frequency,  Ob  is  the  angle  between  B  and  the  x-axis,  and  ay  is  the  polarization 
phase  difference.  By  rotating  the  whole  system  around  the  y-axis  by  the  elevation  angle  yy.  and  then  ro¬ 
tating  that  around  the  x-axis  by  the  azimuthal  angle  f9/c,  we  get  the  full  equation  for  the  magnetic  field  at 
the  center  of  the  cube: 


B  =  B 


cos  Ok  cos  tpk  cos  Ob  cos(cut)  —  sin  Ok  sin  Ob  cos(cut  +  ay ) 
sin  Ok  cos  tfk  cos  Ob  cos(cut)  +  cos  Ok  sin  Ob  cos(wf  +  ay ) 
—  sin  ipk  cos  Ob  cos  (cut) 


The  flux  received  by  each  sensor  should  correspond  with  the  component  of  this  vector.  Parameter  t 
represents  the  time  shift  a  signal  will  undergo  based  on  the  distance  between  a  given  detector  and  the  cen¬ 
ter  of  the  cube: 


c 

where  s  is  the  position  vector  pointing  from  the  cube  center  to  the  sensor,  c  is  the  speed  of  light,  and  k  is 
the  normalized  Poynting  vector.  Since  each  sensor  is  shifted  from  the  center  in  the  same  direction  it  points, 
the  x  component  of  k  is  used  for  the  x  sensor  (in  spherical  coordinates).  Lastly,  a  moving  average  is  used 
to  clean  up  remove  the  higher  frequencies  in  the  output  and  shift  the  data  vertically  so  that  it’s  approxi¬ 
mately  oscillating  around  zero. 

To  simulate  the  electric  field,  equations  for  the  B  field  were  clockwise  by  90  degrees  in  the  polariza¬ 
tion  ellipse,  so  that 


E  = 


cos  Ok  cos  Lpk  sin  Ob  cos  (tut  +  ay)  +  sin  Ok  cos  Ob  cos  (tut  +  ay) 
sin  Ok  cos  (fk  sin  Ob  cos  (tut  +  ay)  —  cos  Ok  cos  Ob  cos  (tut) 

—  sin  <pk  sin  Ob  cos  (tut  +  ay) 


(6) 


Figure  5  generates  Equation  (6).  To  reduce  error,  the  k  values  when  one  of  the  E  or  B  vectors  are  small 
are  removed  before  taking  a  simple  average. 


Poynting  Vector  Calculation  Using  Polarized  Light 


Figure  5.  Both  electric  and  magnetic  fields  rotating  within  the  polarization  ellipse. 


5.2  SIMULATIONS 

Direction  of  arrivals  from  the  Poynting  vector  method  are  very  accurate;  they  have  less  than  a  degree 
of  error  overall,  which  improves  with  longer  detection  time.  The  Poynting  vector  method  also  exhibits 
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good  accuracy  for  varying  frequency  and  polarization.  Figure  6  shows  the  error  in  the  Poynting  vector 
for  a  linearly  polarized  source.  The  simulation  length  was  set  as  three  wavelengths  long  and  the  elevation 
was  fixed  to  45  degrees.  Each  sensor  was  separated  from  the  center  by  one  half  of  a  centimeter  while  the 
highest  frequency  of  3  GHz  corresponds  to  a  wavelength  of  10  cm.  Ideally,  each  of  these  three  sensors 
would  colocate,  but  the  separation  only  has  a  significant  effect  on  accuracy  when  it  has  the  same  order  as 
the  wavelength  or  shorter. 


Figure  6.  SQUID  Poynting  vector  DF  for  a  linearly  polarized  source. 


Using  the  same  simulation  setup  but  a  circularly  polarized  source,  yields  the  results  plotted  in  Fig¬ 
ure  7.  The  circularly  polarized  light  case  gives  a  greater  error,  but  this  effect  is  small  for  longer  wave¬ 
lengths. 


Figure  7.  SQUID  Poynting  vector  DF  for  a  circularly  polarized  source. 
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6.  CONCLUSION 


This  technical  report  offers  an  introduction  into  the  outcomes  of  applying  DF  algorithms  to  data  ob¬ 
tained  from  an  array  of  bi-SQUID  array  sensors.  We  showed  that  it  is  possible  to  perform  DF  tasks  with 
bi-SQUID  arrays,  though  it  is  reasonable  to  assume  the  accuracy  of  the  DF  algorithms  with  a  true  bi-SQUID 
array  will  vary  to  some  degree.  Techniques  presented  here  focus  on  a  design  that  is  as  small  as  possible  to 
maximize  the  number  of  applications  for  which  the  device  can  be  used. 

We  determined  that  the  MUSIC  algorithm  is  a  flexible  subspace  method  that  can  perform  in  various 
situations.  In  addition,  MUSIC  can  perform  DF  with  multiple  signals  of  the  same  frequency;  however, 
MUSIC  requires  k  +  1  sensors  for  k  signals.  If  signal  polarizations  arc  circular  or  elliptic,  MUSIC  can¬ 
not  reliably  find  the  DoA  angles.  The  DR-MUSIC  algorithm  can  determine  the  DoA  of  signals  featuring 
varied  polarization  from  the  bi-SQUID  array  simulated  V (t).  The  DR-MUSIC  algorithm  performed  better 
with  smaller  baselines  than  the  MUSIC  algorithm. 

Excellent  polarization  results  occur  when  a  cube  of  perpendicular  bi-SQUID  arrays  paired  with  a  bi¬ 
directional  electric  field  antenna  using  the  Poynting  vector.  To  ignore  the  separation  between  sensors,  the 
wavelength  must  be  longer  than  a  certain  length,  which  is  outside  the  wavelengths  we  arc  exploring  in  this 
study.  Electrical  field  antennas  may  restrict  the  usage  of  the  DF  device  due  to  their  size. 

Since  bi-SQUID  arrays  can  have  multiple  arrays  on  a  single  chip,  the  projects  team's  goal  is  to  deter¬ 
mine  the  possibility  of  performing  an  accurate  DF  when  arrays  differ  in  location  by  a  matter  of  microns. 
Future  improvements  and  signal  processing  techniques  will  be  implemented  to  boost  the  accuracy  and  util¬ 
ity  of  these  techniques  for  bi-SQUID  array  sensors.  In  this  report,  we  conduct  computer  simulations  of 
bi-SQUID  array  sensors  in  lieu  of  a  real  array  data.  Future  work  will  apply  existing  DF  methods  to  experi¬ 
mental  data  taken  from  a  system  of  SQUID  (or  bi-SQUID)  arrays. 
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